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TEAR-OFF VERSUS GLOBAL EXISTENCE FOR A STRUCTURED MODEL OF 
ADHESION MEDIATED BY TRANSIENT ELASTIC LINKAGES. 

VUK MILISIC AND DIETMAR OELZ 


Abstract. We consider a microscopic model for friction mediated by transient elastic linkages introduced 
in lanQ]. In the present study we prove existence and uniqueness of a solution to the coupled system under 
weaker hypotheses. The theory we present covers the case where the off-rate of linkages is unbounded but 
increasing at most linearly with respect to the mechanical load. 

The time of existence is typically bounded and corresponds to tear-off where the moving binding site 
does not have any bonds with the substrate. However, under additional assumptions on the external force 
we prove global in time existence of a solution that consequently stays attached to the substrate. 
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1. Introduction 

Adhesion forces at the cellular and intra-cellular scales play an important role in several phenomenons 
as cell motility (see [12] and references therein), or cancer growth [14]. In [121 E] the authors derive a 
complete model for a moving network of actin filaments polymerizing near the boundary of the cell and 
depolymerizing close to the nucleus, providing biologically plausible steady-state configurations of the cell 
shape. The main advantage of this method is that the parameters used are experimentally easy to obtain 
if not already available in the literature [a s siiziiis]. The adhesion and the stretching between filaments 
are written as friction terms obtained through a formal limits of a delayed system of equations. Indeed, let 
e be a dimensionless parameter denoting the ratio of the typical lifetime of bonds vs. the overall timescale 
of the model. The asymptotic limit is obtained assuming that the rate of linkage turnover becomes great as 
well as the stiffness of the bonds (typically as 0(l/e)). The rigorous justification of the limit e —?> 0 is the 
ultimate goal of our investigations [ano]. Nevertheless the highly non-linear nature of the delayed model 
leads to consider already the case of a fixed value of e. In this article we show that, even then, the data of the 
problem determins the well-posedness of the model : the balance between the on-rate of the linkages and the 
force exerted on the adhesion point is essential. Mathematically this is seen since, according to this balance, 
either we can show blow up in finite time or global existence. Physically this means that pulling too strong 
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the adhesion site causes a tear-off, and that our model is able to reproduce this feature. Experimentally this 
is observed and used in order to measure the probability distribution of unbinding forces [Hill [mil]. 

More precisely, this study is concerned with a system of equations which describe the evolution of the 
time-dependent position of a single binding site as it moves on a ID-substrate. External forces / act on the 
moving point-object while it is attached to the substrate through continuously remodeling elastic linkages 
which represent the effect of transiently attaching protein bonds. Their age distribution is denoted by 
Pe = Pe{t,a) where a > 0 denotes the age of linkages and t > 0 denotes time. Here we treat e as a fixed 
constant, which we keep in our notation to maintain consistency with previous studies, in which we were 
analyzing the convergence with respect to s 

In [121 111] the following structured model for the turnover of protein bonds with age distribution p^ = 
Pe{t, a) for age a > 0 has been established, 

{ edtPe + daPe + C Pe = 0, t > 0, a > 0 , 

Ps(t,a = 0) = p,o) , t>0, 

Pe{t = Q,a) = ppe{a), a>0, 

where po{t) := pg{t, a) da and the on-rate of bonds is a given coefficient j5 times a factor, that takes 
into account saturation of the moving binding site with linkages. This system is coupled to the elongation 
Ug = Us{t,a) of the linkage through 

(1.2) C := C(|we(t,a)|) , 


according to which the off-rate ^{u) > 0 is a real, positive function of the elongation of the linkage Ug. In 
m we introduced the following age-structured model for the evolution of the elongation 

edtUe + daUe = — ( edtf -I- [ {C{ue)uePe) {t, d) dd] , t > 0 , a > 0 , 
do \ Jo J 

Ue{t, 0 = 0) = 0, t>0, 

Ue{t = 0, o) = Ui^e{a) , O > 0 , 

where / = /(<) S ffi. is the external force acting on the binding site. In [10] we had shown that the system 
(11.31) is equivalent to an integral equation for the position of the binding site Zg = Ze{t) itself |9], namely 



(1.4) 



Ze{t) 


{Ze{t) 

-- Zp{t) 


- Z^it - ea)) Pe{t, a) da = f{t) , 


t > 0 , 
t < 0 , 


where the known past positions are given by the Lipschitz function Zp(t) G R for t < 0. The correspondance 
between dm and dm is made through 

Zsit) - Ze(t - sa) 

Us — 

£ 

Note that the equation m has been the original result of the mathematical modeling since it represents 
the balance of external and elatic forces acting on the binding site. On the other hand, in [10] it has turned 
out to be beneficial to work on the system (11.111 . (11.311 instead, since it allowed to derive powerful a-priori 
estimates on Ug. 

The analysis in the older studies [5] and nni relied on the existence of an upper bound Cmax of the function 
C- In [13, for fixed e, we prove existence and uniqueness of weak solutions of the coupled system dm, m, 
m- In the simpler semi-coupled case when (^(t, a) is a given bounded function we give as well a convergence 
result as e —?> 0. The analysis of dm in [3 relied on the fact that ^ 'Ho[p —Poo] < —Cmin/£’Ho[p —Poo] where 


(1.5) 

JiolpKt) := [ \p{t,a)\da + 

/ p{t^a)da 


Jr+ 



and where poo is the stationary solution of dm- For the analysis of (11.311 in [10] we established the a-priori 
estimate (Lemma 4.1, [lOj l 


ldtf(t)ldt 


( 1 . 6 ) 


Ps(t,a)lue(t,a)lda < 


pi,e(a)lu/r(a)lda + 
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Note that both results, decay of the functional (11.51) and the a-priori estimate (11.61) . do not rely on the 
existence of an upper bound of C, and therefore do hold in the framework of this paper. 


It is the aim of the present study to relax the hypothesis of boundedness of This represents a major 
improvement of the analysis, because the lower bound of the total mass po{t) strongly depends on Cmax and 
the analytical arguments in [10] do rely heavily on this control. Furthermore the upper bound Cmax had 
major importance in the fixed point argument used in m to prove the global existence result since we used 
it to control the non-linear right hand side in (11.31) . 

In addition to deepening the analysis, unboundedness of the off-rate is the natural scenario from the 
application point of view. A typical situation is Bell’s law, i.e. an exponential increase of the off-rate as the 
elastic linker is extended, ( = Coexp(|u£|) (cf. [THIS])- However, this strongly non-linear scenario is still 
out of reach of the rigorous mathematical analysis that we present in this study which relies on C being a 
(globally) Lipschitz continous function as it’s main technical assumption. 

The right hand side of dESl for a given function u, 

9u{t) \edtf[ C{u{t,a)) gu{t,a)u{t,a)da 

do,u ( Jb+ 

where Qu solves (HH) with ( = C(u) and /io,« := Quit, a) da, can become infinite if either ^o,u vanishes or 
/r+ ^ du da blows up. We define the modified right hand side 

= _ 1 

■ max(^o,ix,M) 

where p, and p are two strictly positive arbitrary constants. The strategy to prove our existence result is first 
to establish existence and uniqueness of a solution of this modified problem using a fixed point argument in 
the space 


edtf max —p, min p, 


Ciu) Qu 




(1.7) 


Xt := 


eLi“((o,r)x 


-) s.t. 


sup ||u(t, a) a;(a)||j;^oo < oo 
te(o,T) 


defined for any specific time T > 0, where the weight function is 


( 1 . 8 ) 


uj(a) 


1 

1 -(- u 


To this end we introduce the map $ : u G Xt i—>■ u G Xt where, given v, we solve dnj with C 
obtain the age distribution . Then we look for the solution of the problem : 


{ edtu-{-dau = ^^{t) , t>0, a>0, 

u{t, 0) = 0 , t > 0 , 

u(0,a) = uj^g(a) , a > 0 , 


((v) and 


to obtain u G Xt- The right hand side of (11.91) becomes a bounded function whose bounds depend on the 
cut-offs p and p. This allows to prove contraction of the map $ on a time interval that is sufficiently small. 
Due to the uniform bounds this process can be iterated to obtain (p^,, w) a unique, global in time, solution. 
Then we establish a uniform bound on p(t) := (^iw)wpu,da, the second integral term in g^,,. This shows 
that for p sufficiently large with respect to 1/p, p{t) never reaches p so that the solution {p.u,,w) satisfies also 
a simple-cut-of problem where can be replaced by defined as : 


9u-= -7-- C{u)Quuda \ . 

ma.x{po^u,P) [ J 

In a second step, we prove that if additional assumptions hold, this solution never reaches the cut-off 
value p. Otherwise, we give a lower bound to the time span during which the cut-off is not reached. In both 
cases the solution of the modified problem is also the unique solution to the original system dnj-dESl) either 
globally in time or on the finite interval of time. 
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More precisely, in Section 21 we analyze the dependence of the lower bound of with respect to the 
L°°(0,T) norm of This naturally leads to local existence results for the original problem (ll.ll) - (ll.3l) in 
Section [S] by providing a minimal time for which the solution (p^,, w) does not reach the cut-off value fi. 

Even stronger results are rigorously obtained in Sections [6] and [7] generalising a straightforward compu¬ 
tation in the special case where ^(u) = 1 -|- |m| and assuming that remains strictly positive. In this case, 
integrating (HD in age, and using the fact that (ED transforms in a)ue{t, a)da = f{t), we obtain 

that 

edtpo - /3(1 - Mo) + Mo + / = 0 , 

which can be solved directly. This provides immediately the bounds 

min (po(0), < Mo(t) < Mo(0) (l " > 

where 

min “1“ 1 V /min - /3 max 

and leads to a strictly positive lower bound of po when /3niin > /max, whereas if fnun > jd-max, the time to 
is well defined and the binding site tears off, i.e. po{t) becomes zero, at t = to- These basic ideas provide 
global existence results (Section [SD vs. tear-off results (Section [7]) under more general assumptions on Q. 


2. Technical assumptions, preliminary results and a-priori estimates 

2.1. Hypotheses. 

Assumption 2.1. a) There exists a minimal value Cmin s.t. C.{u) > Cmin > 0, Vw G K. 

b) The derivative of C, is bounded i.e. |C^(w)l ^ Clip, Vu G M. 

c) The function f is Lipschitz continuous on [0,T] for any positive fixed T. 

Remark 2.2. Note that with this definition we do not allow more than a linear growth for But in contrast 
to [HI [in] one has no hypothesis concerning boundedness on (f. 

As in m we assume also some hypotheses on the initial and boundary data of (inD: 

Assumption 2.3. The initial condition is 

(i) nonnegative, i.e. pi.eifl) > 0 , a.e. m R+ . 

(a) Moreover, the total initial population satisfies 

nOC 

0 < / ppsia)da < 1 , 

Jo 

(Hi) and higher moments are bounded, 

poo 

0 < / a^pi^e{a) da < Cp , for p = 1, 2 , 

Jo 

where Cp are positive constants depending only on p. 

Assumption 2.4. For (3 we assume that 

a) P = P(t) is a continuous function, 

b) 0 < Plain < Pit) < Pmax for all positivc timcs t. 

We detail hereafter those results from [S] that are still valid in the weaker frame of Assumptions 12.1112.31 
and 12.41 


Theorem 2.5. We suppose that u is a given function in Xt. Let Assumvtions \2. li \2.3\ and \2.4\ hold, then 
for every fixed s there exists a unigue solution g G C'°(IR+; L^(]R+)) fl L°°(]R^) of the problem (11.11) . with the 
off-rate Cit,a) := C(R(t,a)). It satisfies (II.ip in the sense of characteristics, namely 


P{t — ea) (l — g{d, t — ea) da) x 

g(t^a) = { X exp (— C(d, t — e(a — d)) dd) , whena<t/e 

Pi,eia-t/s)eyip(^-J J*Ciii-t)/sa,i) dtj , ifa>t/e. 


( 2 . 1 ) 
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Lemma 2.6. Under the same assumptions as in Theorem \2.5\. let g be the unique solution of problem (11.11) . 
then it satisfies a weak formulation of this problem, namely 

nOO pT poo 

(2.2) / / g{t,a) {sdt^pda^ — C^) dtda — e / = T,a) da-\- 

Jo Jo Jo 

pT poo 

+ / g{t,a = 0) ip{t,0) dt + e / pi^e{a)ip{t = 0,a) da = 0 , 

Jo Jo 

for every T > 0 and every test function (p G T] x R-i-). 

Following the same argumentation as Lemma 2.2 in one has 

Lemma 2.7. Under the same assumptions as in Theorem \2.5[ it holds that po{t) < 1 for any time. This in 
turn implies that g{t, a) > 0 for almost every {t, a) in 

For p G N we define the p-th moment of the solution of (HH) 

poo 

Ppit) := / a^g{t,a) da . 

Jo 

Then, following the same argumentation as Lemma 2.2 in [^, one has 
Lemma 2.8. Under the same assumptions as in Theorem \2 . 51 

h’pij) ^ /^p,max for p = 1, 2 , 

where the generic constants /ip,max read : 


P'- /3n 


/ip,max:-^ . I G 

i^o ^'Cmin Pmm + (,n 

Proof. When p = 0 we simply integrate (ED with respect to age : 

edtpo + Ppo + / fpeda = jS 
Jr+ 

as C is bounded from below and using Gronwall’s lemma one has : 

/-^max 


Uoit) < po(0) 


+ Cn 


-'min I Smm 


For any integer p we then write : 


^bJtpp “b Cmin/ip PTp—1 — b 
P 


which using Gronwall’s lemma gives 

\\Tp\\l’=°{o,t) < /ip(0) + ||/ip_l||p,oo^p 

Smin 

By recursion, one proves the claim. 


□ 


Proposition 2.9. Under Assumvtions \2.1[ 1 2., 91 and 12.41 setting p ■= g 2 — gi where p 2 and gi solve ED 
with off-rates C(w 2 ) (resp. ((wi)) where W 2 (resp. wi) is a function in Xt, we find that 

'Ho[p]{t) < co(l - exp(Cmint/e))||u'||jf^ , Vt G (0,r) , 

where w := W 2 — wi, cq := t-^ClipA^i max; /ii.max being the bound on the first moment of gi. 

Proof. The proof follows the same lines as for Lemma 3.2 and Lemma 3.3 in [5] based on the system satisfied 
by p, 

edtp + daP + C 2 P = -Cqi t > 0, a > 0, 
p{t,0) = —f3{t) / p{t,a)da, t > 0, 

J 

/3(0, a) = 0, a > 0, 


where ( := ((^ 2 ) — C(wi). 


□ 
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For /c > 1 we define 


'Hk[p]'= [ {1 + a)''p{t,a)da 
J M-i_ 


for these functionals one has : 

Proposition 2.10. Under the same hypotheses as in the previous proposition, and if moreover 


( {1 + aYpi^e(.o) da < oo, V£g{ 0, fe + 1}, 

M-L 


then 


'Hk[p]it) < hfc(l - exp(-CminV^))||w||^^, Vt G (0,r), 

where the constants hk depend only on Cmin)CLip and on the constants (/i^,max)^G{o.fc+i} related to the bound 
on the i-th moment of Q 2 - 

Proof. We apply a recursion argument. The case fc = 0 is proved by Proposition 12.91 We suppose that the 
claim is true for £ < /c — 1. We have formally that 

£dt{l + a)^|p| + 9a(l + a)^|p| ~ ^(1 + ^\p\ + Cmin(l + a)^|p| < |CI(1 + aY §2 

Integrating in age, one gets that 

edinklp] - /3|/t| + Cmin'Hfc[yo] < knk-iip] + CLip||w|lxt / (1 + aY+^g2{t,a)da 

Jr+ 

which is then estimated giving: 

edt'Hk[p] + Cmin'Hkip] < kHk-llp] + CLipC'fc+l||w|lxt + /3maxHo[p] 
which using the Gronwall’s Lemma gives 
1 6Xp( fminl/ 


Hfc[/5](t) < 


Cn 


sup {kUk-ilpKs) + /3max'Ho[p](s) + CnipCfc+i ||u;||^ ) 
e(o,t) 


where we used, in the last estimates, the recursion hypothesis and Proposition [2i9 


□ 


If we give ourselves T > 0 and a function g G L°°(0,T) and then we compute w as the solution in the 
sense of characteristics of 


(2.3) 


{ edtw + daW = g(t) , f > 0, a > 0 , 
w{t, 0) = 0, t > 0 , 

r(;(0, a) = ui^^{a) , a > 0 . 


And all along the paper we will assume that the initial condition uj^e belongs to L°“(K+,w). For this simple 
transport problem it holds that 

Theorem 2.11. If T > 0 and g is a function in L°°{0,T), for any fixed e and any T > 0 there exists a 
unique w G Xt solving problem (IQ) . Moreover one has the a priori estimates: 

Moreover the maximal time of existence is infinite if g € L°° (R+). 

3. Global existence results for cut-off problems 
We solve the problem find {q,w) satisfying : 

' edtQ -I- daQ + C,{uj)pe = 0 , t > 0, a > 0 , 


(3.1) 


Q{t, 0) = /3(t) j a)da^ , 


^ £>(0,0) = pi,e{a) , 


t > 0 , 
a > 0 
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and 

(3.2) 

where we set 

(3.3) f^(t) := 


{ edtw + daW , t>0 

= 0 , 

w(0,a) = u/,e(a) 


a > 0 , 
t>0 , 
a > 0 , 


1 


fi) 


edtf + max —p, min 


C{w)pewda, p 


where /io(t) = /g Pe(t, a) da. The two constants p, and p are positive. 


Lemma 3.1. We suppose that {p,p) £ 

GiA,B) := 


5.^)^ and that f is Lipschitz. The function 
{edtf + max(-p, min(p, B))} 


max(/i, A) 

is a Lipshitz function with respect to £ R for any fixed i? £ R and with respeet to S £ R for any fixed 
A £ R. The Lipschitz constants in both cases are uniform and depend only on {p,p). 

Theorem 3.2. We suppose that Assumptions \2.1l and 


hold. Moreover we assume that £ 
L°°(R+,a;) and is finite and that the constants p andp are fixed. For any fixed time T possibly 

infinite, there exists a unique pair of solutions {gw,w) £ (7(0, T; L^(R_|_)) x Xt solving the coupled problems 
(CT . (lO) and (lO) . 

Proof. We apply the Banach fixed point Theorem to mapping w £ Xt i-t m £ Xt such that 

£dtu +daU = g,„{t), t>0,a>0, 

w(t, 0) = 0, t > 0, 

w(0, a) = ui^e(a), a > 0, 

We prove that $ is actually contractive in Xt for a time T small enough, 
a) The map $ is endomorphic. For any given w £ Xt one has invariably 

(3-4) WJ < - (4dtf\\i 

p V 


lL“(0,T) +Pj J 

which by the same method as in Theorem 12.111 provides a bound independent on T in Xt on u : 
b) The map $ is a contraction. We set g^ ■= gw 2 ~ 9wi p := — Qwi and so on. Thanks to Lemma 


=:/i+/2. 


\9w(t)\ < |£||5t/|lL«>(0,T) +p} + ~ 


'1 1 

f A 

f+ - 

/ Qgwda 

J p 

\Je+ / 


Ii is immediately estimated thanks to Proposition 12.91 and one has : 


h < 


^ {e|l'9t/llL-(0.T) +p}'^o[p]it) < {e||3t/||ioo(0,T) +p}co||w|U,, 


while we decompose the difference of triple products in I 2 as : 
1 


I 2 < — 

a 

1 

< - 


CQw2W2 + C 1 PW 2 + CiQwiWda 


CLip|R’|Pu;2|a^2Ma 


(CLipIkillxt+Co) < / {l + af\p\da\\w 2 \\x,+ j 4 + af Qwtda\\w\\x^ 

I J ]R_i- J K_i_ 


< cl hi 


AT. 


B.2[p\{f)) < c||w||^^, 
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where the constant c depends on Clipj Co) (l|w^i||js:j)ie{i, 2 }j and a^pi^e(a)da for k G {0,1, 2} . Using 
again Theorem 12.Ill one has 


lx. < 


t 


t 


tc, 


II5u)|Iloo(Q < ll5l«|lLoo(0,t) — ll^^llxt' 


t + e' e" "" £ 

If To < e/c then there exists a unique fixed point w G Xtq of the mapping $. 
c) Global existence for any time. We suppose that existence is established on the whole time interval 
[0, Tn-i] for n > 1. We construct a fixed point for the next interval [T„_i, T„ := T„_i + AT„] on the map 
u = $(?;) 

( edtu + dau = t G (T„_i,T„),a > 0, 


and 


<1 u(t,0) = 0, 

[ u(Tn-i,a) = w{Tn-i,a) 

' £dtp + daP + C{v)p = 0, 

p{t,0) = P{t) ^1-J p{t,a)da^ 


t G (T„_i,T„), 
a > 0. 

t G (T„_i,T„),a > 0, 
t G (T„_i, T„), 
a > 0. 


Pe{t,a) := 


We ■■ = 


w{t,a) if t e [T„_i,T„) 
w{t,a) tG{0,Tn-i] 


^ piTn—i) a) — piTn—i, a), 

If we denote the extensions to [0,r„] of (p,u) as : 

if < e [T„_i,T„) 
g{t,a) tG{0,Tn-i] 

where w = ^{w) and {g,w) is the unique solution of (I,1.II) - (I3.2|) on [0,T„_i]. The continuity of pe allows 
to apply Lemma l2.8l Similarly for Wg one has 

\\We\\xT„ - I|5uj(0X[T„_i,T„) +9w^X[Q,T^-i]\\xT^ + IIl~(R) 

^ (£||at/llL~(0.T„) +P). 

^ - - - + lk/,e|lLj,(R). 

tL 

where XA is the characteristic function of the set A, and we used the uniform estimate on provided 
by (13.41) . These estimates prove that the constant c in the contraction in b) is not changing as time 
evolves. Thus we can fix-point again choosing AT„ as in the previous paragraph and prove contraction 
in [T„_i,T„]. At this step the recursion is complete. The theorem is proven for any positive time. 

□ 

Corollary 3.2.1. Under the same hypotheses as above, for any pair of positive definite reals {p,p), the 
solution-pair (g^^jw) solving (IXTD-(IO) satisfies the a priori estimates (fm . 


Proof. We use that 


\9wit)\ < 


1 


< 


Uo{t) 

1 

Mo(t) 


£\dtf\ min p, 


({w)wgw da,p 


£\dtf\ -f min J Ciw)gn,\w\da^ | < |e|9t/| -f J Ciw)gni\w\do 


Then the same arguments as in the proof of Lemma 5.1 in m provide the a priori estimates. Indeed in the 
sense of characteristics |w| satisfies : 


£9*1^1 -f 5a|w;| <\9w\< 


1 


poit) 


e| 5 t/l + / Ciw)g.u,\w\da } . 


then multiplying the later inequality by and integrating with respect to age, one gets : 


edt 


gw\w\da 


C{w)\w\gnjda < £\dtf\-G / C{w)\w\gwda 


and because on the right and on the left hand sides the same integral terms cancel, the claim follows. □ 
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Proposition 3.3. Under Assumvtions \2. R \2.3\ and \2.4\ let {qw,w) be the solution of the fully coupled and 
stabilized problem (IO)-(IO)-(l331). there exists a positive finite constant 72 s.t. 

/ <^{w{t,a))\w{t,a)\gw{t,a)da < —, Vt > 0 , 

where the constant 72 is depends on 

• the a priori bound only on gi^\w\da (obtained in Corollarv \3.2.1\) , 

• Clip, and C(0). 

Proof. Using equations (EH), EH and hypotheses 12.11 one has 

edt{gw\w\C) + 5a(£'u,|?u|C) + < gw\w\{edtC + 9oC) + CQw\Vw\ ■ 

Integrating in age and setting p(t) := gwit, a)\w{t, a)\f{w(t, a)) da gives 


edtp+ C'^\w(t,a)\gw(t,a)da < \gj i Chip gn,\w\da+ Ciw)gw{t,a)da\ 

< \Vw\ (^Cup da + C(0)^ 

< -{£\dtf\+p) (2CLip/7i+C(0)), 
tL 

where £iu,|w|da < I/ 71 . Now we consider the second term in the left hand side above: using Jensen’s 
inequality one writes 

f lR^Ciw)\w{t,a)\gw(t,a)da\‘^ ^ J^^{Ciw))‘^\w{t,a)\gu,{t,a)da 


\ \w\g^da 

since |?u|pu,/ /r |w|£»uj da is a unit measure. This implies that 


Jr \w\gyjda 


{Ciw))‘^\w{t, a)\g^{t, a) da > 


(/r+ C{w)\w{t,a)\gyj{t,a)da'j 


> HP ■ 


/r_^ \w\g^,da 

We obtain a Riccati inequality 

edtp + HP^ <h/tL + p/li, P(0) = / C('W/.e(a))|u/,e(a)|p/,e(a) da , 


where h := e||9t/||^ (2CLip/7i + C(0)) is a constant. We denote by P± the solutions of the steady state 
equation associated to the last inequality, i.e. P solves jiP'^ — P/g. — h/p, = 0. The solutions are given by 

1 


-P± = - (1 ± + 4 / 1 ^ 71 ^ /( 271 ) < - max (|p( 0 ), (^1 ± y/l + 4 / 171 ^ /( 271 )) =: 


72 
M ' 


Applying Lemma lA.il we conclude that p{t) < max{p(0), P-p} < h/Ui which ends the proof. 


□ 


Theorem 3.4. Suppose that Assumvtions \2. 1[ lil.dl and \2.4\ hold, moreover, suppose that € L°“(]R+,a;) 
and that ll^t/llioo^g t) Is finite, if {pw^w) is the unique solution of the stabilized problem (I3.1D - (I3.2D - (I3.3I) . it 
is also the unique solution of EH-EH together with the modified right hand side : 

(3.5) 5 ^ =- ^ - -(edtf+ [ C{w)wgruda] . 

ma.x{go,w,g) y Jr^^ J 

Proof. The proof is a simple application of the Proposition 13.31 above and taking p > 72 /^ when solving 
(I3.1ll - (j3.2l) - (j3.3|) . Indeed, in this case one has that the truncated right hand side from (13. 3 1) becomes (13.5L 
one since p(t) := C{w)wgyjda never reaches ±p. □ 
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4. Impact of the cut-off value on the mean bonds’ population 

In this section we give ourselves a function g S L°°(0,T) and compute w S Xt solving (12.31) . In what 
follows we analyze the properties of an age structured model for linkages whose off-rates depend on w : we 
define g^j as the solution of 

{ edtQw + daQw + ((w) Qw = 0, t> 0 , a >0 , 

g^(t,a = 0) =/3(t) (1 - go,w) , t> 0 , 

g^{t = 0, a) = pi,e{a), a > 0 , 


where go,w{t) := gw{t,a)da. 

We compute a sharper upper bound on namely 

Lemma 4.1. Let Assumvtions \2. 1\ and \2.4\ hold. Let g^ he the solution of (14.11) . We suppose that /io,uj(0) < 
1. Let us fix a positive constant 70 s.t. 

70 < min ("1 - /io.^(O), -—^ . 

\ Smin “r Pmax / 

Under Assumvtions \2. 1[ \2.A and \2.4\ Uo,w{t) <1 — 70 holds for every positive time t. 

Proof. We proceed similarly as in Lemma 2.2 in [9]. The computations are thus only formal although they 
can be made rigorous exactly as therein. By hypothesis, the data satisfies 1 —70 —Aio,iu(0) > 0. By continuity 
this also holds on a time interval [0,to) small enough. We proceed by contradiction and suppose that at time 
to the mass po,w(to) reaches 1 — 70 . The equation on /io,u, reads: 


edtpo.w - - Uo,w) + / gw{t,a)C,{w{t,a))da = 0 . 

Multiplying it by —1 and estimating /3(t) < /3niax, one deduces that 


(1 70 UO . w ') 4” /5max(l 70 PO . w ) 4” 7o/^max 

then the lower bound on implies 



gwit,a)C.{w(t,a))da > 0 , 


70 UO.w) 4“ /3max(l bO M0,u>) 4- 7o/^niax ^ CniinM0,ia ■ 

We transform the latter right hand side writing 

CminMo.u’ ~ Cmin{l TO UO.w) Cmin(l bo) ■ 

Setting q(t) := (1 — 70 — go.wft)), one then has 

4- (Cmin 4- ft max )<?> Cniin (Cmin “1“ ^max )7o > 0 , 

the latter estimate being true under the hypothesis that 70 < Cmin/(Cmin 4- /3max)- The conclusion then 
follows integrating the latter inequality in time 

/"j- \ ^ ( (/3max 4- Cminj^O \ ^ 

q[to) > exp I -^- I g(0) > 0, 

under the hypothesis that 70 < (1 — fJ.o,w{0))- But this contradicts the assumption that q{to) = 0, which 
ends the proof. □ 


We do not have a positive definite lower bound on /io,u, yet : at this stage we only know that po^wft) > 0. 
For this reason we define gf^(t,a) := g(t,a)/{iu.o,w{t) + S) and we observe that this new function is in 
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L“^((0,T) X R-I-). It solves the equation 


(4.2) 


edt'Qi, + dagi, + 


-P 


giit,a = 0) = Pit) 


1 


Mo,iu 

1^0,w ^ 1^0,w “t“ ^ 

1 f^0,w 


Mo.ui + S i-io,w + / 


^giit = 0, a) = pi,eia)/iPo,w + 5 ), 


g„ = 0, t > 0, a > 0 , 

t>0 , 
a > 0 . 


When (5 = 0 one denotes gf^ simply by g^. 

Proposition 4.2. Let g € L°°{0,T) be given, and let (g,w) he the solutions of (14.II) - (12.31) . Under Assump- 
tions \2.l\ 13.31 and \2.4\ and if pLo,wiO) < 1~70; there exists a constant^ independent of S and e such that for 
every positive 6 it holds that 


tR+ 


C(w(t,a)) giit,a)da < C + Cu 


lip lly iIl°o(o,t) 


loP, 


T 

mm ^ 


Vt > 0, 


where C, := C(0) + C,iui,ei^))Pejio)da. Taking the limit as d goes to 0, one obtains then the analogous 
result for . 

Proof. The product p{t, a) := Qg% satisfies 

edtp + daP + ^j = Ciw)9it)gi ■ 

Indeed, using arguments as in Lemma 2.1 p. 489 and Lemma 3.1 p. 493 [5], one proves that if w solves (ESD 
and f is Lipschitz, then ({w) solves {edt + da)C{w) = C'{w)g in the sense of characteristics (as in Theorem 
12.51) with the corresponding boundary conditions. Then the latter equation on p is understood in the same 
manner. 

Integrating in age and setting g(t) := pit, a)da we conclude that 

(4-3) £dtq-Cit,0)gl,it,0)+ [ Cgida-([ + 9el(i, 0) < CLipllsIloo • 

yt/K_|_ J 

To find a lower bound for p^(0,t) we choose 6 < 7o/2 and use the upper bound on po,wit) established in 
Lemma Q] in order to obtain 

1 


(4.4) 


giiOU) > Pn 


_ 11 > « . 22 . 

^ I ^ Mmm 2 • 


,1-70 + 5 

Assuming po^wit) > 0 we also find using Jensen’s inequality that 

([ Ciwit,a))giit,a)da\ <( iCiwit,a)))^g^da .x < [ iCiwit,a)))^gida . 

\Jr+ J Jr+ iho.w + d) 

If Po,wit) = 0 the same inequality holds true since then git, a) = 0 for almost every a. These considerations 
allow then to rewrite (14.31) as 

£dtq + gii0,t)iq - CiO)) < CLipllfflloo • 

Setting q := q — C(0) and using Gronwall’s Lemma gives 


qit) < exp ( -i / giiO, s)ds'] g(0) + 


exp - 


1 


Thanks to the uniform lower bound (14.41) we conclude 


qit) < exp - 


P min'lpt 
2e 


m 


2CLip||g|L 

70/dmin 


1 — exp — 


J gii0,s)ds] dr . 
Pminlot'' 
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which then gives turning to the variable q that 
(4.5) qit)<C{0)+[ C{ui^sia))pej{a)da + 


2CLip||ff|l(x) 


1 — exp ( — 


le 


./K+ loP min 

This bound is uniform in S. One passes to the limit 5 = 0 which gives the final result. 

Proposition 4.3. Under the same assumptions as above and if po,w(t) < 1 — 7o o-nd choosing /lo.min s.t. 

Pn 


□ 


/^0,min < min /io.u)(0), 


+ C + Cui] 


lin ( —J 

y 'ToPmin ^ ^ J . 


NyilL-(0,T) 

where we used the bound provided by Proposition \4-S\ one has a lower hound on ■ 

PQ,w{t) > Po ,min 1 vt > 0 . 

Proof. We integrate the equation (ED with respect to age which gives 


£dtPQ,w - P{i - Pq,w) + gw{t,a)C{w{t,a))da = 0 , t>0 

JR+ 

Po,w{0) = / pi,e{a)da , 

J ]R_L 


t = 0 . 


we divide and we multiply the last term on the left hand side by po,w and we write: 

edtpo.w + ^/ 3(0 + J CPwda^ po,w = P{t) , t > 0 . 

Now suppose that there exists a time small enough s.t. po,w{t) > Mo,min for all t € [0, to) and that po,w(to) = 
Mo,min- We use the notation A := ^ + CLip I|5|Il“(o t) i^ -yoP ’ ?) write for the difference Mo(f) := 
Mo ,W (t) - Mo ,min 


“t“ I Pif) P I CPwda j Mo — /^(0(f Mo, min) I / CPwda j Mo, min ^ 

y Jr+ j \Jr+ j 


^ /5min(l Mo, min) '^M0,min — (/^min “t“ 


Pn 


Mo, min I 7 0 , 


^ Praia + A 

which holds thanks to the bound on da established in Proposition 14.21 and the definition of mo,i 

Using Gronwall’s Lemma we finally obtain that 

PO.wifo) Mo,min 7 exp f f {^P{t) A)iiT | (MO,nj(0) Mo, min) 7 0 , 


which contradicts the fact that po,w{to) = Mo, min- This ends the proof. 


□ 


5. Local existence of the fully coupled problem 

Theorem 5.1. Let f be a Lipschitz function on (0,r) and G L°°(K+,a;). We suppose that Assumptions 
roiro and \2.4\ hold. Let {gw,w) be the solution of ED-(ED together with the simple cut-off defined 
by (13.51) . Then for any fixed p < po,w{Q) there exists a time 

T = - {PraiaP ~ {Praia + C)M^) 

73 - 

for which po,w{t) > M for any t G (0,T). So the solution {g,w) of (I3.1I) - (I3.2I) - (I3.5I) is also the unique local 
solution of the fully coupled system ED-iED- 

Proof. Gathering results above one has : 
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since we suppose that ^ < 1 and we set 73 := CLip(e||9t/||^oo(o t ) + 72 )- Thanks to ProDOsition l4.31 the lower 
bound on then becomes : 

/ N . f /3minAI 

> mm I fj,o,w{0), — - 


(/3min + 


Choosing fj, < /ro,u)(0) we tune T s.t. 




(/3min + 


73'T 


> M 


□ 


6. Global existence eor specific data 


Under hypotheses of Theorem 13.21 whatever be the time of existence T for the solutions of the 

stabilized model, then thanks to Corollarv l3. 2.11 one has that : 

/ C{w{t,a))Q{t,a)da < / (C(0) + CLip|w|)pu,da 


< C(0) + Clip 


\ui^e\pi,eda + / \dtf\ds] =: (, Vt G (0,T). 


Proposition 6.1. Under assumvtions \2.1[ and \K^ if (3rain > C set : 

c 


0 < Aio.min < min 1 - 


MQ,w(0) 


then 


has 


Uo,wit) > Mo .mini vt e (o,T). 

Proof. We set ft := pLo,w{t) — Mo,min and write the equation that it satisfies : 

Pp> — / fpda /3(1 Mo,min) ^ C ~b /5min(l M0,min)- 

•/M-f 

We estimate from below the right hand side using previous bounds. The lower bound is positive definite 
provided that Praia > C and that Mo,min < 1 — C,/Praia- Using Gronwall’s Lemma, one has : 

p{t) > exp(-/3maxi/e)M(0) > 0 

if Mo,min < Mo,u;(0), which ends the proof. □ 

Theorem 6.2. If we fix a finite time T > 0. Under Assumvtions 1X7] and \2.tA and assuming that 
i) f is Lipschitz on (0,T), 

ii) P satisfies assumptions \2.4\ together with /3min > C 

there exists a unique solution {pg,Ue) € (7(0,T;L^(K.+)) x Xt solving system (ll.l|) - (ll.3l) . 

Proof. By Theorem 13.41 there exists a unique couple (pu,,ic) G G(0, 00 ; L^(]R+ )) X Xoo solving (I3.1I) - (I3.2I) - 
(1331) for any given constant p. We choose T > 0 and provided that P satisfies hypothesis required by 
Proposition lB.ll we set the constants 0 < m < Mo,min according to Propositions 16.II Then mo,™ does not reach 
the threshold value p so that 


9w{t) = 


1 


max(Mo,-i«(i),M) 
1 


sdtf + / {C{w)gww){t,a)da \ = 

k R4- / 


/^0,ip (0 


sdtfP {Ciw)gn,w)(t,a)da\ = gm{t), a.e.tG{0,T). 

J 


The pair {gyj, w) is in fact also solving (ll.ll) - (ll.3l) on this time interval. This provides existence of a solution 
{pe,Ug) = {gw,w) on [0,r]. Since by Theorem 13.41 (o,„. w) is unique, so is {pe,Ue) in this time period. □ 
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7. Blow up for positive solutions 

Theorem 7.1. Under assumvtions \2.S\ and if Tq is the time of existence of{pg,Ug) solving (|1.1I) - (I1.3I) . and 

i) ui^e{a) > 0 for a.e. a € IR.+ , 
a) dt 'f{t) > 0 for a.e. t e (OjTq), 

then the product pg(t,a)ue{t,a) is non-negative for a.e. {t,a) € (0,To) x IR+. 

Proof. Since it holds that /(O) = pi^s(a)ui^g{a) da and f{t) = pg(t,a)ue{t,a) da yields 


Peit,a)\ueit,a)\da < / p 7 ,£(a)|M 7 ,£(a)|(ia + / \dtf{t)\di = 


/ ppe{a)ui^e{a)da-\- / dtf{i)di=f{t)= / pe{t,a)ue{t,a) da, 
JBj- Jo JRi 


which implies the result. 


□ 


Lemma 7.2. Suppose that fc is a convex locally differentiable function. Then for any function u G X^o, 
one has : 

C'(0) / upeda< / Cciu{t,a))pe{t,a)da - CciO)po, a.e. t G M+ 

Proof. Since Cc is convex, for almost every {t, a), one has : 

C'(0)(?7(t, a) - 0) < Cciu{t, a)) - C(0) 

and integrating with respect to Peda, one has the desired result. □ 

Proposition 7.3. Under assumvtions \2.‘J\ and \2.4\ and if 
i) C satisfies Assumvtions 1X71 and admits a lower convex envelop fc s.t. fc(u) < Ci^) for all u G K+ with 

CM > 0, 

a) let f be a Lipschitz function s.t. dtf{t) > 0 for a.e. t G (0,T), 

Hi) f and fi are s.t. < CMfmin, 
iv) ui^e{a) > 0 for a.e. a G K.+ , 

then if the solution {pe,Us) solving (ll.ll) - (ll.3l) exists until a finite time Tq, this time cannot be greater than 


to := 


In 1 + 


Mo(0)(/^min + Cc(0)) 
CMM - /3max 


/3min + Cc(0) 

for which 

Tofto) < 0 . 

Moreover, on (0,to) x R+, one has a lower bound on the profile of Ug namely 

, . , min(t,ea) 

Ug{t,a) > £76 In I 1 + 


{to-t) J’ 

where ye := to infte(o.to) dtf/poiO)- 

Proof. By Theorem l7.1l n^ft. al > 0 a.e. {t,a) G (0,To) x R+. The equation for po reads : 

edtpo - 13{1 - po )/ ({ue{t,a))pe{t,a)da = 0 
that we estimate using Lemma 17^ as follows : 

edtpo - PC - Po) P CA'J) I Ue{t,a)pe{t,a)da + Cc{0)po < 0 
and becomes under these simplifications : 

(7.1) sdtPo — /3(1 — po) + (c(O)po + CcOf ^ 0- 

We can deduce from this equation that 

edtpo + (/3min + Cc(0))po P P max Cc(0)/min 
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which gives using Gronwall’s Lemma that fJ-oit) < M(t), where 
p(t) := /ro(0) exp + Cc(0)) _ C'MUin - /3„ 


1 _ exp ( _ 


(/3min + Cc(0)) 

Looking for the time to s.t. 7l(^o) = 0 provides the explicit form of to in the claim. Thus To < to. Moreover, 
as 7l(t) is a convex function one has that ; 


t \ — 


to 


t 


Mo(i) < 1 - — m(0) + —Kto) = 1 - :r Ai(0) 


to 


t \ 


to 


and because, by Lemma [7.21 C,{ue)uePe is positive almost everywhere on (0,to) x R+, 

sdtf ^ £76 


cdtU^ “t“ daUe ^ 


> 


a.e in (0,to) x 


Po(t) to -1 

Using Duhamel’s formula provides 

^ if t > EU, 

^Ui,e(a - t/e) + £76 otherwize, 

which then gives the lower estimate on u^. 




□ 


Appendix A. Riccati inequalities 

Lemma A.l. Let £ > 0 and real, let y be a positive differentiable function oftG K+, satisfying 

f edty + Ay^ <By + C, t > 0, 

\ 2/(0) =2/0, t = 0 

where yo > 0 and {A,B,C) £ (M-i-)^. Setting y+ := {B + ^B'^ + 4AC)/(2A), one has that 

y(t) < max(j/o,2/+), Vt G R+. 

Proof. We set m := max(yo, y+), it satisfies —Amf + Bm + C < 0. Then we define y := y — m which then 
solve the differential inequality : 

(A.l) edty + Ay"^ + {2mA-B)y <0, 

Since the quadratic term is positive we neglect it, and apply Gronwall’s Lemma : 


y{t) < exp - 


{2Am — B)t 


y{Q) = exp - 


{2Am — B)t 


{yo - m) < 0 


which ends the proof. 


□ 
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